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ABSTRACT 

We present an entirely analytic model for a preheated, polytropic intergalactic medium in hydro- 
static equilibrium within a NFW dark halo potential in which the evolution of the halo structure 
between major merger events proceeds inside-out by accretion. This model is used to explain, within 
a standard ACDM cosmogony, the observed X-ray properties of nearby relaxed, non-cooling flow 
groups and clusters of galaxies. We find that our preferred solution to the equilibrium equations 
produces scaling relations in excellent agreement with observations, while simultaneously accounting 
for the typical structural characteristics of the distribution of the diffuse baryons. In the class of 
preheating models, ours stands out because it offers a unified description of the intrahalo medium 
for galaxy systems with total masses above ~2 x 10 13 M©, does not produce baryonic configurations 
with large isentropic cores, and reproduces faithfully the observed behavior of the gas entropy at large 
radii. All this is achieved with a moderate level of energy injection of about half a keV, which can 
be easily accommodated within the limits of the total energy released by the most commonly invoked 
feedback mechanisms, as well as with a polytropic index of 1.2, consistent with both many observa- 
tional determinations and predictions from high-resolution gas-dynamical simulations of non-cooling 
flow clusters. More interestingly, our scheme offers a physical motivation for the adoption of this 
specific value of the polytropic index, as it is the one that best ensures the conservation after halo 
virialization of the balance between the total specific energies of the gas and dark matter components 
for the full range of masses investigated. 

Subject headings: cosmology: theory — galaxies: clusters: general — methods: analytical — inter- 
galactic medium — X-rays: galaxies: clusters 



1. INTRODUCTION 

In this paper we explore the consequences of apply- 
ing a recently developed procedure for the evolution of 
the halo structure onto the modeling of the properties 
of the hot gaseous baryons they host. Our model as- 
sumes, in line with the results of numerical simulations 
(see Salvador-Sole, Manrique, & Solanes 2005), that be- 
tween major mergers halos grow inside-out maintaining 
invariant their structural parameters. This concept has 
been the basis for the successful reproduction of both the 
NFW (Navarro, Frenk, & White 1997) functional form 
of the dark matter density distribution and the empirical 
correlations involving its concentration parameter (Man- 
rique et al. 2003; Salvador-Sole et al. 2005), while it has 
also been invoked to resolve some apparent challenges to 
the CDM model for structure formation (Loeb & Peebles 
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2003). We are interested in confronting the predictions 
of this halo growth process with both the observed global 
and structural X-ray properties of nearby groups and 
clusters of galaxies and, in particular, in investigating 
its compatibility with the findings of both observations 
and high-resolution gas-dynamical simulations that the 
diffuse baryons outside the innermost regions of galaxy 
systems obey a polytropic equation of state that is nei- 
ther purely isothermal nor adiabatic. 

It has been known for nearly two decades that X- 
ray galaxy systems are at variance with the simple self- 
similar scenario predicted by Kaiser (1986) in which grav- 
itational collapse is the sole mechanism driving the evo- 
lution of their nonbaryonic and baryonic components. 
Provided that galaxy groups and clusters of total virial 
mass M can be treated as scaled versions of each other 
and that their X-ray emitting gas is in hydrostatic equi- 
librium within the dark matter (DM) potential, this see- 
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nario assumes that the virial temperature 2 of the dark 
halo, T* ir = GM/2r v ir, is a good approximation for the 
effective X-ray gas temperature, T x = kTx/p> = <r|, 
where r v i r is the virial radius, ax the ID velocity dis- 
persion of the gas, and p, — /im p the mean molecular 
weight of the gas in physical units, /j, being the cor- 
responding value in units of the proton mass m p and 
k the Boltzmann's constant. Under these conditions it 
readily follows that T x oc M 2/3 and, if the total X-ray 

2 /3 

gas mass M x oc M, that T x oc M x ' too. Further- 
more, if the plasma emissivity is dominated by thermal 
bremsstrahlung — this is strictly true only for the hotter, 
non-relativistic plasmas with Tx > 3 keV ~ 3 x 10 7 K, 
as at lower temperatures, emission line cooling by metals 
dominates — , its integral over the entire energy range of 
the (continuum) X-ray emission and over the gas distri- 
bution, assumed to be an isothermal plasma with (num- 
ber) density nx = Px/ft proportional to the mass density 
distribution and extending likewise out to r v ; r , predicts 
a total halo bolometric X-ray luminosity Lx oc T x or, 
equivalently, Lx oc M 4 / 3 . 

While hydrodynamic simulations dealing only with 
physical processes that do not have any preferred scale 
(e.g., Navarro, Frcnk, & White 1995; Evrard, Metzler, 
& Navarro 1996; Bryan & Norman 1998; Eke, Navarro, 
& Frenk 1998; Frcnk et al. 1999; Thomas et al. 2001) 
have extensively confirmed the approximate validity of 
the self-similar scaling even for cosmogonies of the CDM 
type, years of observational efforts have now provided 
compelling enough evidence for an excess in the spe- 
cific energy/entropy of the X-ray gas with respect to 
the contributions arising from adiabatic compression and 
shock heating alone (e.g., White 1991; David, Forman, & 
Jones 1991; David, Jones, & Forman 1996; Ponman, Can- 
non, & Navarro 1999; Lloyd-Davies, Ponman, & Cannon 
2000; Finoguenov et al. 2002; Ponman, Sanderson, & 
Finogucnov 2003). These deviations are indicative of an 
important role of nongravitational processes in determin- 
ing the properties of the hot intrahalo medium (hereafter 
IHM). 

Basically, there are two different approaches competing 
to explain the observed X-ray properties of bound collec- 
tions of galaxies and, in particular, the deviations from 
the simple self-similar picture. One category of models 
(both analytical and numerical) explores the effects of 
nongravitational heating on the IHM drawing inspiration 
from the work of Kaiser (1991), who was the first to point 
out that the intracluster gas properties, including the 
evolution of its X-ray luminosity function, could be bet- 
ter reconciled with the hierarchical scenario of structure 
formation if the entropy of the gas we now see was the re- 
sult of some early injection of heat — perhaps during the 
epoch of galaxy formation — previous to the cluster as- 
sembly. The addition of energy into the gaseous compo- 
nent increases its entropy and reduces the shock-heating 
efficiency, preventing the gas from reaching high density 
during the halo collapse and producing lower than ex- 
pected luminosities for a given temperature (e.g., Cava- 
liere, Menci, & Tozzi 1998; Wu, Fabian, & Nulsen 1998; 
Balogh, Babul, & Patton 1999; Locwcnstcin 2000; Wu, 

2 In the present work, temperatures expressed in specific energy 
units will be identified by a star superscript. 



Fabian, & Nulsen 2000; Tozzi & Norman 2001; Bialek, 
Evrard, & Mohr 2001; Borgani et al. 2001; Babul et al. 
2002; McCarthy, Babul, & Balogh 2002; Dos Santos & 
Dore 2002; Oh & Benson 2003). For a given amount 
of energy AE injected into the intergalactic medium, 
this extra heating is expected to leave a more signifi- 
cant imprint on the low-temperature systems for which 
AE/E > 1 and the shock-created entropy is negligible, in 
agreement with observations. This means that for these 
systems gas infall proceeds roughly adiabatically, pre- 
dicting the existence of an 'entropy floor' (Lloyd-Davies 
et al. 2000). Paradoxically, radiative cooling is also ca- 
pable of explaining the lack of self-similarity of the X-ray 
gas (e.g., Bryan 2000; Voit & Bryan 2001; Muanwong et 
al. 2001; Wu & Xue 2002; Dave, Katz, & Weinberg 2002). 
In this alternative scenario, the central low entropy gas, 
characterized by a cooling time shorter than the typ- 
ical halo age, is selectively removed by condensing into 
dense, cold structures and rapidly replaced by the higher- 
entropy material from the halo outskirts. The latter is 
then heated by adiabatic compression as it flows in, giv- 
ing rise to temperature profiles with a central maximum 
and engendering substantially shallower density distribu- 
tions that cause a reduction in the bulk X-ray luminosity, 
much like in the non-gravitational heating scheme. 

These two basic scenarios, however, are not free of 
important drawbacks. Radiative cooling without non- 
gravitational feedback mechanisms is a runaway pro- 
cess that results in severe overcooling of gas on group 
scales and thereby in an unpleasantly large fraction of 
gas (<~ 50%) converted into a 'stellar' cold medium (Oh 
& Benson 2003). This so-called 'cooling crisis' or 'cool- 
ing catastrophe' (e.g., Blanchard, Valls-Gabaud, & Ma- 
mon 1992; Fukugita, Hogan, & Peebles 1998; Balogh 
et al. 2001) is difficult to be dealt with appropriately 
even with state-of-the-art numerical simulations due to 
the technical challenge represented by the simultaneous 
requirement of properly implementing small-scale pro- 
cesses, such as star formation and feedback, and having 
a large dynamic range from galaxies to large-scale struc- 
tures in a volume large enough to guarantee a statisti- 
cally representative ensemble of simulated halos. Over- 
coming overcooling is not easy either within analytical 
radiative cooling schemes (e.g., Voit et al. 2002). Im- 
portant cooling and condensation of the gas, and the 
internal heating resulting from the ensuing star forma- 
tion feedback, continuously modify the gas distribution, 
so the necessary hypothesis that the IHM is in hydro- 
static equilibrium does not hold. Besides, any heating 
of the intergalactic medium after it is confined within 
halos of groups and clusters has the disadvantage of re- 
quiring a larger amount of energy injection to raise the 
denser gas to a given entropy level (Loewenstein 2000; 
Voit & Bryan 2001). Precisely, the major criticism to- 
ward the preheating scenario has to do with the 'energy 
crisis' (e.g., Dave et al. 2002) related to the yet unidenti- 
fied astrophysical source(s) responsible for injecting the 
amount of excess energy (<~ 1 keV particle -1 ) necessary 
to explain the observed scaling relations and its connec- 
tion with the process of galaxy formation (e.g., Navarro 
et al. 1995; Cavaliere et al. 1998; Balogh et al. 1999; 
Lloyd-Davies et al. 2000; Wu et al. 2000; Brighenti & 
Mathews 2001; Babul et al. 2002). We want to point 
out, however, that although important details on how 
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and when this additional heat is injected into the diffuse 
medium have still to be sorted out, there is no shortage 
of potential sources of pre-collapse heating (see, for in- 
stance, Babul et al. 2002 for a brief summary). More 
importantly, the substantial additional entropy at large 
radii detected in both clusters and groups of galaxies 
at a typical level of <~ 400keVcm 2 (Finoguenov et al. 
2002; Ponman, Sanderson, & Finoguenov 2003) appears 
to be beyond the bounds of plausibility of models that 
invoke only gas cooling and gravitational heating to ex- 
plain the ~ 100 hr 1 ^ keVcm 2 entropy excess seen in the 
central regions of galaxy systems (Ponman et al. 1999; 
Lloyd-Davies et al. 2000; Ponman et al. 2003). Two re- 
cent results give convincing support to this conclusion. 
To begin with, there are the hydrodynamic simulations 
by Finoguenov et al. (2003) with radiative cooling, star 
formation, and non-gravitational heating, which demon- 
strate that star formation without extra heating pro- 
duces too steep entropy profiles in group-sized halos that 
also fall somewhat short in explaining the height of the 
observed specific entropy on the halo outskirts. In addi- 
tion, observations of the hot cluster Abell 1795 by Ikcbc, 
Bohringer, & Kitayama (2004) show that the intraclus- 
ter medium temperature is higher than the dark matter 
'temperature' (measured from the radial profile of the 
dark matter velocity dispersion), even in the central re- 
gion where the radiative cooling time is short. 

While it is true that recent high spatial resolution X- 
ray data suggest that treatments combining energy injec- 
tion with cooling may be necessary to explain the cen- 
tral temperature and entropy profiles of massive cooling 
flow clusters (e.g., McCarthy et al. 2004), this work aims 
to demonstrate that a hydrostatic polytropic gas model 
based on the notion of preheating and requiring a rela- 
tively modest heating energy budget can have a remark- 
able success in accounting simultaneously for the global 
and structural X-ray properties representative of the bulk 
of present-day relaxed galaxy systems with gravitational 
masses between ~ 10 13 and ~ 10 16 M©. This is done 
in the context of the concordance flat ACDM cosmol- 
ogy with reduced Hubble constant h — 2/3, present-day 
matter density f2 m = 1/3, rms mass fluctuation on scales 
of 8 h~ x Mpc a% — 0.95, cosmological baryon content 
Ob = 0.04 (and, hence, a cosmic baryon fraction of 0.12), 
and primordial mass density fluctuation power spectrum 
index n = 1, fully consistent with recent joint analy- 
ses of WMAP and redshift survey data (e.g., Tegmark 
et al. 2004). The outline of the paper is as follows. In 
Section 2, we provide a brief overview of the approach 
used to follow the growth of bound halos and present 
the profiles describing the structure and kinematics of 
their nonbaryonic and baryonic components. After iden- 
tifying the best solution to the equilibrium equations of 
the hot baryons in Section 3, we then go on in Section 
4 to validate our model predictions for nearby (z = 0) 
galaxy systems against an extensive set of X-ray observa- 
tions involving a variety of both scaling relations between 
bulk properties and structural properties of the IHM, 
discussing the results in the light of recent observations. 
The last section of the paper contains a summary of our 
main findings. 



2. THEORETICAL FRAMEWORK 

We present here the set of equations that completely 
defines our model for the evolution of the structure of 
bound X-ray groups and galaxy clusters. As usual in 
analytic approaches of this kind, these systems will be 
approximated by spherically symmetric structures whose 
global dynamics is exclusively driven by the assemblage 
of the dark matter component (baryons amount < 15% 
of the total gravitational mass M) . 

2.1. The Dark Halo Evolution Model 

Numerical simulations (e.g., Navarro et al. 1997; Jing 
2000; Bullock et al. 2001) of the hierarchical mass as- 
sembly in the universe, have demonstrated that the mass 
density distribution in the ACDM cosmology of relaxed, 
nonbaryonic halos at any redshift z can be characterized 
by a universal NFW functional form 

r 3 

P(r) = Pc , ° , (1) 
r(r s + r) 2 

where p c and r s are, respectively, the halo characteristic 
density and scale radius. 

Salvador-Sole, Solanes, & Manrique (1998; see also 
Raig, Gonzalez- Casado, & Salvador-Sole 1998) developed 
a consistent analytic description of the hierarchical evo- 
lution of DM halos that provides an excellent fit to N- 
body simulations (Raig, Gonzalez-Casado, & Salvador- 
Sole 2001). This framework for structure formation in- 
corporates in the well-known extended Press-Schechter 
formalism (Lacey & Cole 1993) a pre-established phe- 
nomenological threshold A m setting the fractional mass 
increase that separates the two basic mass aggregation 
regimes of dark halos: minor mergers or accretion, where 
these objects grow inside-out through the continuous ag- 
gregation of small clumps that do not disturb their in- 
ternal structure 3 , and major merger events, in which the 
progenitor halos are fully disrupted giving rise to the 
formation of a new bound system. While, as shown by 
Raig et al. (2001), it is possible to achieve an excellent 
agreement between theory and simulations for any value 
of A m , it is essential specifying the exact value of this 
parameter in order to explain the characteristic inner 
structure of halos (Manrique et al. 2003), as well as any 
correlation involving their concentration (Salvador-Sole 
et al. 2005). In particular, the concentration dependence 
on halo mass at any redshift found by Eke, Navarro, & 
Steinmetz (2001), which we shall adopt here, requires 
A m = 0.21 for the elected cosmology (Hiotelis 2003). 
Once A m has been chosen, the epoch of formation of a 
population of halos of mass M — defined as the redshift 
at which they have experienced their last major merger — 
is fixed. We take the median value of the analytical prob- 
ability distribution of formation times for DM halos with 
observed virial masses M (see eqs. [6]-[12] in Salvador- 
Sole et al. 1998) as their typical formation redshift Zf OT . 

By defining x — r/r s as the radial distance in units of 
the scale radius, we can rewrite equation (1) in dimen- 
sionlcss form as 



3 This condition implies that the values of p c and r s remain 
unaltered during the inside-out growth process. 
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where, taking into account the relationship 



M = —r 3 vil A virPu 



(3) 



with A v j r the ratio of the mean density of a sphere of 
radius r v ; r to the characteristic background density of the 
universe p u at the epoch in which halos are observed, one 
has that p c = A v i r p u c 3 /3[ln(l + c) — c/(l + c)], with c = 
r v ir/r s the halo concentration giving the reduced virial 
radius of the system within which the mean halo density 
is A v ; r /9 U . In our model, p u = p CI it(z) = 3H 2 (z)/8irG 
and A vir (z) = 178f^ 45 (z), resulting in A vir w 100 at 
z = 0. 

This dimensionless form of the halo density profile is 
used to find the dimensionless mass within radius x 



~, . M(x) , .„ . x 
M{x) = -^y 1 = ln(l + x) - 



(4) 



l + x 1 

the corresponding dimensionless halo circular velocity 
square 



V\x) = 



V 2 {x) _ ln(l + a;) 1 



V 2 



1 



(5) 



and, from the Poisson equation r s d$(x)/dx = GM(x)/x 2 
and expression (4), the dimensionless NFW halo poten- 
tial 

$(x)^=-^i±^, (6) 

where the (invariant after halo formation) character- 
istic parameters M c = 4nr 3 p c , V 2 = GM c /r s , and 
$ c = —$(0) = V 2 act as normalization constants. 

Besides, for the computation of the halo total energy 
(§ 3.1) it will be also necessary to deal with both the 
reduced one-dimensional velocity dispersion of the dark 
matter, which we obtain numerically by integrating the 
Jeans equation, assuming isotropic orbits and a null pres- 
sure at infinity, 



a 2 {x) EE 



aHx) 
V? 



1 



p(x) 



p(x) M ^ dx . 



(7) 



and the masses halos with present-day total masses M 
had at Z{ OT , which we infer from the rate of halo mass 
growth by accretion given in Raig et al. (2001). 

2.2. Modeling the Hot Gas Component 

2.2.1. Dimensionless Profiles 

The equilibrium structure of the IHM is determined 
by requiring this component to be, in epochs of gen- 
tle mass accretion, in thermal pressure-supported hydro- 
static equilibrium within the (fixed) gravitational poten- 
tial wells set by the dark matter halos 4 . Hydrodynamical 
simulations on the physics of diffuse baryons in accretion 
flows indicate that most of the bulk energy of the in- 
falling gas is converted through highly efficient shocks 
into thermal energy and that the outer boundary of the 
hot gas expands following closely the growth of the col- 
lisionless dark matter component (see, e.g., Takizawa & 

4 An implicit and necessary hypothesis of this kind of models 
is that both the gas and the dark matter swiftly readjust to a 
new hydrostatic equilibrium after a major merger event; see, for 
instance, the discussion in Cavaliere, Menci, & Tozzi (1999). 



Mincshige 1998). So, as a practical approximation for 
the evolution of the IHM structure, we will assume that 
between major mergers the hot gas distribution evolves 
inside-out maintaining its outermost radius permanently 
equal to r v j r . Furthermore, we will consider that all the 
baryonic matter within the virialized halos is in the hot 
X-ray-emitting phase (observational data indicate that, 
typically, only < 10% of the baryons in clusters are locked 
into stars and cold gas so fi g ~ f^; see, e.g., Lin, Mohr, 
& Stanford 2003; Bell et al. 2003). 

The Jeans equation for an ideal gas, Px — nxkTx = 
pxTx, with isotropic pressure takes the form 

d*(x) 



d\n Pg (x) dlaT*(x) 



dx 



(8) 



dx dx 
Thus, to solve the IHM structure by relating univocally 
the gas density and temperature to the halo mass more 
information in the form of a relation between T* and p g 
is required. Analytic models obtain conveniently flexible 
solutions by adopting for the IHM profiles a polytropic 
equation of state P cx p 1 (see also § 3.1), leading to the 
relationship 

p^x)^Tl(xfl^\ (9) 

where 7 is the polytropic index — which for a given 
potential effectively specifies the shape of the temper- 
ature profile — that is presumed independent of the halo 
mass and radial distance to the halo center. Since we 
are adopting the common approximation of describing 
the hot X-ray plasma by a single-phase (i.e., each vol- 
ume element contains gas at just a single temperature), 
monoatomic perfect gas with a null metallicity gradient, 
for a strictly adiabatic IHM 7 must be equal to 5/3, the 
ratio of specific heats. Gas models with a polytropic 
index larger than this value are convectively unstable, 
while isothermal distributions are retrieved for 7 = 1 
(apart from the cores of some clusters, there is little ob- 
servational evidence for radially increasing X-ray tem- 
peratures, corresponding to 7 < 1). 

The solution of the equilibrium equation (8) by ap- 
plying boundary conditions at the current value of the 
virial radius of the halo (in r s units), i.e., at x = c, 
leads to the following dimensionless temperature profile 
for < x < c : 

f g ^)ESM = f g *(c) + I-l($( c )-$(,)) . (10) 

Therefore, from equations (9) and (10) the dimension- 
less gas density writes 



Pg(s) 



s (*) 



Pg(c) 



p g (c)exp 



1 



In 



Tj(x) 

7-1 

7 



1/(7-1) 



In 



>(*) 



where p g 
gas and 



Pg,iso(c) 



,(H) 



p c f2b/^m is the characteristic density of the 
$(c) - $(ar) A 



Pg.isoOc) = Pg,iso(c) exp 



T*(c) 



(12) 



represents the solution of the Jeans equation for an 
isothermal gas with a dimensionless temperature whose 
value is equal to T*(c). Equations (10)-(12) completely 
determine the structure of the hot gas once its tempera- 
ture and density at x = c are specified. 
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2.2.2. Boundary Conditions 

To determine the value of the additive constant T*(c) 
in equation (10), we have taken into account the fact 
that, in the absence of non-conservative processes, the 
total specific energies of the gas and dark matter must 
remain invariant and equal to each other, as numerical 
simulations confirm (Thomas ct al. 2001; Muanwong et 
al. 2002) — we are neglecting the small gas heating that 
may arise from any energy transfer between the dark 
matter and the hot gas occurring during gravitational 
collapse (see, e.g., Pierce, Thomas, & Couchman 1994). 
Thus, when the preheating energy of the gas is non-zero, 
we have, for any redshift z < Zf OI : 

£ g = £ DM + AE* . (13) 

In equation (13) £i = (K + U)i/Mi represents the ra- 
tio between the sum of the total kinetic and potential 
energies and the total mass of the i component cal- 
culated within the reduced virial radius c(z), while in 
AE* = AE/p, pre the non-starred AE is chosen to rep- 
resent from now on the excess energy per gas particle 
brought by preheating and p pre the mean mass of the 
gas particles where this energy is injected. 

By substituting in equation (13) the model profiles de- 
rived in the current and former sections, we obtain, after 
some algebra, the implicit relationship 

J C (3T*(x) + $(*)) p g (x)x 2 dx - 2AE*/V 2 
j^p g (x)x 2 dx 

^(?>a 2 {x)+$>{x)}p{x)x 2 dx 

J^p(x)x 2 dx ' ^ 

which we choose to solve at the observed redshift z = 
(see also § 3.1). In this last equation the value of ^ pre is 
chosen by assuming that energy injection takes place af- 
ter reionization is complete (z < 6), but before the galaxy 
systems have formed (recent estimates by Finoguenov et 
al. 2003 and Oh & Benson 2003 suggest that preheat- 
ing may have happened at z < 3, a value substantially 
larger than the typical formation redshifts < 1 our model 
predicts for group and cluster-sized halos). Accordingly, 
we consider that the excess energy is injected into a fully 
ionized plasma with metallicity equal to the current stan- 
dard 30% solar abundance and temperature above 10 5 K 
— so p becomes temperature independent — for which 
Mprc = 0.998 x 10~ 24 g, as inferred from the tables of 
Sutherland & Dopita (1993). Note that, for the smallest 
systems, the nongravitational heating energy AE is ex- 
pected to be the best part of the energy budget, making 
their mean gas temperature almost independent of the 
system mass. 

To set the normalization p g (c) of the dimcnsionlcss gas 
density profile (eq. [11]), we choose to impose the bound- 
ary condition, also adopted in other models with pre- 
heating (see, e.g., Oh & Benson 2003), p g (r v i r ) / p(r v i T ) = 
^b/^m implying that 

p g (c) = p(c) . (15) 

The results of both adiabatic simulations and simula- 
tions with cooling and star formation (Ascasibar et al. 
2003; Kravtsov et al. 2005) show that the local ratio of 
gas to dark matter density within present-epoch halos is 



close to the cosmic baryon fraction near the halo bound- 
ary. Therefore, by using equation (15), we implicitly 
assume that, in the range of halo masses where our pre- 
heating model is aplicable (see § 3.2), the baryon frac- 
tion near the virial radius is not strongly affected by the 
non-gravitational heating. Besides, a 'local' constraint 
— as opposed to the more widely used 'global' bound- 
ary condition M g = (Qb/Q m )M — allows for the pos- 
sibility that the finite central temperature of the IHM 
(sec, e.g., cq. [10]), which causes the central gas distribu- 
tion to be more flattened than the density distribution 
of the dark matter, can produce baryon fractions within 
the virial radius smaller than the cosmic value. Indeed, 
given that any energy injection into the gas actually en- 
hances this tendency by creating a core in its density 
profile and decreasing further its central value, equation 
(15) grants that our model can reproduce the reduced 
total gas masses characteristic of low-mass halos in pre- 
heating scenarios (Muanwong et al. 2002; Oh & Benson 
2003; see also the last paragraph of § 4.1). 

3. BEST SOLUTION TO THE EQUILIBRIUM EQUATIONS 

Before we investigate the X-ray properties of real 
galaxy systems with our model, it is necessary to identify 
the correct solution to the equilibrium equations of the 
hot baryons. In our analytical treatment, the gas struc- 
ture is controlled by the values of the polytropic index 
7 and the preheating energy per particle AE, which we 
now proceed to determine. The appropriateness of the 
resulting best model will then be exhaustively checked 
in Section 4 by confronting its predictions with an ample 
variety of X-ray data on both the global properties and 
radial structure of the gas. 

3.1. Fixing 7 from the Specific Energy Balance 

While 7 appears to be a parameter relatively well con- 
strained empirically, the consensus on its value is by no 
means absolute. A number of observations (e.g., Marke- 
vitch et al. 1998; Finoguenov, Arnaud, & David 2001; 
Finoguenov, Reiprich, & Bohringer 2001; Sanderson ct 
al. 2003) and both analytical models (e.g., Cavaliere et 
al. 1999; Tozzi & Norman 2001) and high-resolution gas- 
dynamical simulations (e.g., Ascasibar et al. 2003; Bor- 
gani et al. 2004; Ettori et al. 2004), agree in showing 
that the IHM of galaxy clusters is well represented, out- 
side any cooling region, by a polytropic equation of state 
with a universal index of around 1.2 ± 0.1 over the full 
radial range. There are, however, some recent observa- 
tions suggesting both that the polytropic treatment of 
the gas is only acceptable beyond about 20% of the clus- 
ter virial radius and that the required value of 7 is close 
to 1.5 (De Grandi & Molendi 2002). As we next show, 
the physical consistency of the inside-out growth of the 
halo structure in regard to the energy balance between 
its two main components can be used to infer a well de- 
fined value for this index, thereby implying that 7 is not 
a true degree of freedom in our model. 

As stated in Section 2.2.2, relation (13), which was 
used to set the boundary condition on the equilibrium 
dimensionless temperature profile, must hold from Zf or 
provided a negligible fraction of the hot gas is able to 
cool after the halo formation. While we have chosen to 
fix the balance of total specific energies between the two 
main halo constituents, £ g — £dm, at the observed redshift 
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Fig. 1. — Model percent deviation at Zf OI between the sum 
of the kinetic and potential specific energies of the hot gas and 
its corresponding total specific energy inferred from equation (13) 
plotted against the total virialized mass of present-day halos. 
Comparison is made for three different levels of energy injection: 
AE = OkeV particle -1 (panel a), AE = 0.5 keV particle -1 (panel 
b), and AE = 1.0 keV particle -1 (panel c). In all panels, different 
line types denote different values of the polytropic index: 7 = 1 
(solid), 7 = 1.1 (dotted), 7 = 1.2 (short-dashed), 7 = 1.3 (long- 
dashed), and 7 = 1.4 (dot-short-dashed). 



z = 0, obviously there is no guarantee whatsoever that 
its universal value, AE*, is preserved automatically for 
any other cosmic time larger than the virialization epoch. 
In fact, as the equilibrium equations of the inside-out 
evolving DM and polytropic gas profiles are different, 
the degree to which condition (13) is fulfilled along the 
halo lifetime depends on the specific values adopted for 
the parameters 7 and AE controlling the gas structure. 

Figure 1 summarizes the (maximum) percent differ- 
ence between the sum of the kinetic and potential spe- 
cific energies of the hot gas and the total specific en- 
ergy of this component inferred from equation (13) cal- 
culated at Zfor- Results for three different levels of mass- 
independent AE, each one for five different values of 7, 
are presented. We note that the physical consistency of 
our predictions has been also preserved by limiting the 
calculations to halo masses for which the sign of the sum 
of the kinetic and potential energies of the IHM is nega- 
tive, as otherwise it cannot remain gravitationally bound 
(see also § 3.2). Figure 1 shows that, in general, the 
deviations increase with decreasing present-day gravita- 
tional mass reflecting the more important contribution 
of the extra heating to the gas energy budget in small 
halos. In panel (a), where the results corresponding to 
a null preheating (AE = 0) are displayed, the value of 
7=1.1 gives the best performance, keeping the relative 
deviation below 1 per cent across the entire range of halo 
masses that are being investigated. The curve for 7 = 1.2 



also shows moderately low deviations ranging from 2% 
at the high mass-end up to somewhat less than 6% for 
the smallest systems. Compared to these two models, 
the performances of the 7 = 1.3 and 7 = 1.4 solutions 
are much less satisfactory. This latter conclusion is also 
applicable to the results obtained for an isothermal gas 
profile (7 = 1) in the low half of the mass range, which 
only leads to results as good as in the 7 = 1.1-1.2 cases 
for halo masses above ~ 10 15 M Q . 

The effects of including nongravitational gas heating 
are illustrated in panel (b) of Figure 1 , were we display 
the deviations calculated for AE = 0.5 keV particle -1 , 
and panel (c), were a larger amount of energy injection 
of 1.0 keV particle -1 has been assumed. The addition of 
preheating energy results in larger discrepancies, except 
for the case 7 = 1.2, which now preserves the universality 
in the balance of specific energies to within 3% across the 
three-order-of-magnitude range of gravitational masses. 
Now, even the behavior of the 7 = 1.1 curves, which is 
highly satisfactory for cluster-scale systems, swiftly de- 
teriorates once they move below ~ 1-2 x 10 14 M©, giving 
rise to unacceptably high energy differences. Under these 
same conditions, the relative deviation of the isothermal 
solutions remains satisfactorily low only for the rarest 
most massive clusters. We note that all curves in Fig- 
ure 1 have been calculated by taking z for equal to the 
median of the formation time distribution of the dark 
halos of a given current mass M. In general, the larger 
the adopted formation redshift, the larger the maximum 
relative deviation that results at this redshift. Thus, for 
instance, by setting Zf OI equal to the lower quartile the 
discrepancies become typically a factor two larger. 

These results all together imply that for a moderately 
preheated IHM 7 ~ 1.2 warrants the best global physical 
consistency with the inside-out growth of structure by ac- 
cretion. Interestingly, polytropic gas distributions with 
7 = 1.1-1.2 in a NFW potential well also have the de- 
sirable property of tracing the dark matter profile in the 
outer part of halos (Komatsu & Seljak 2001). Besides, 
as we discuss in the next section, a universal polytropic 
index of 1.2 for the diffuse baryons provides the best fit 
of our predictions to the observations when X-ray data 
are used to adjust 7 and AE simultaneously. 

3.2. Constraining AE from the 
Luminosity- Temperature Relationship 

The X-ray luminosity and temperature are the most 
easily observable bulk properties of the hot gas in clus- 
ters and groups of galaxies. Accordingly, they have been 
measured for a large number of these systems. This large 
dynamic range of the data and the few hypothesis in- 
volved in the calculation of these two observables, which 
are reflected in a relatively tight correlation, render this 
relationship the most reliable empirical reference for the 
calibration of any theoretical modeling of the IHM. In 
this section, we fix the poorly constrained value of AE 
by finding the theoretical solution that, for 7 = 1.2, best 
matches the observed luminosity-temperature relation- 
ship for groups and clusters of galaxies. 

Comparison to observations is done via the total bolo- 
metric X-ray luminosity, derived by integrating the vol- 
ume emissivity of the X-ray-emitting gas, e(x), out to 
the halo maximum radius: 
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Fig. 2. — Luminosity-temperature relation versus observations. 
Different datasets are identified with different symbols. The solid 
line shows our best-fit model prediction drawn for a polytropic 
gas with 7 = 1.2 and mctallicity equal to 0.3 Zq, which yields 
AE = 0.55 keV particle -1 . A Hubble parameter of h = 2/3 has 
been applied to both the model and the data. 
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and the emission-averaged X-ray temperature of the halo 

expressed in energy units 



V c 2 J^jl(T*,Z)T^x)e(x)x 2 dx 
Jq e(x)x 2 dx 



(17) 



In equations (16) and (17), p s (x) and T* (x) are, respec- 
tively, the mass density and temperature profiles of the 
IHM given by equations (10)-(12), whereas A(T,Z) is 
the bolomctric cooling function normalized to the total 
number density of all plasma species (electrons and ions) 
taken from the tables of Sutherland & Dopita (1993), 
which include both bremsstrahlung and cooling from 
metal lines, for a plasma of metallicity Z = 0.3 Zq. In 
these calculations, we have kept track of the cooling time 
at c, 

3 2 P g (c)f g *(c) 
<-cool — 2 e(c) ' 

in order to exclude halos for which all of the gas within 
the virial radius has had time to cool since the system 
formed. 

Figure 2 presents the comparison of the observed 
Lx — Tx relation for groups and clusters of galaxies 
with our best-fit model, which is obtained for AT — 



(18) 



0.55 keV particle -1 . In this figure, the theoretical curve 
encompasses halos of present-day virial masses ranging 
from 6 x 10 15 M large enough to explain the hottest 
system data, down to 2 x 10 13 M Q . The exclusion of sys- 
tems with masses below the latter value is due to the 
impossibility for halos this small of retaining their IHM 
when energy injection reaches 0.55 keV particle -1 . Inter- 
estingly the emission-weighted X-ray luminosity associ- 
ated with this lower-mass limit (<~ 10 41 ergs s -1 ) shows 
a good correspondence with the minimum luminosity of 
the measurements in the GEMS galaxy group project by 
Osmond & Ponman (2004) when this dataset is restricted 
to those systems with genuine group emission (G-samplc; 
open pentagons in the figure). 

Figure 2 also illustrates that the slope of our best-fit 
model increases gradually with decreasing Tx matching 
the trend of the observations at cluster and group scales. 
In the hot cluster regime (Tx > 3 keV) Lx oc T^ 3 , 
as suggested by most authors (e.g., White, Jones, & 
Forman 1997; Arnaud & Evrard 1999), although our 
model indicates that the logarithmic slope should ap- 
proach asymptotically the self-similar value of 2 for the 
very hottest systems. For the coolest systems, the typ- 
ical model behavior (Lx oc T£ 5 ) is likewise fully con- 
sistent with observations, including the latest data from 
the GEMS project G-sample that, nonetheless, show a 
noticeably increase in the scatter about the predicted 
trend for Lx < 10 42 ergss -1 . As argued by Osmond & 
Ponman, this marked raising of the dispersion — arising, 
for the most part, from pushing the current instrumental 
capabilities to the limit in order to observe a sufficiently 
representative number of cool groups (Tx < 0.7 keV) — 
introduces several sources of bias that likely conspire to 
flatten the Lx — Tx relationship at these scales and ren- 
der it to apparently follow the cluster trend, instead of 
the characteristic logarithmic slope of 4-5 found for X- 
ray bright groups (Helsdon & Ponman 2000; Xue & Wu 
2002). Still, the fact that our model goes through the 
middle of the locus occupied by the data points in the 
low Tx, low Lx regime suggests that the Lx— Tk relation- 
ship does indeed steepen substantially in these extreme 
scales. 

It is also interesting to note that had we decided not 
to fix for a start the value of 7 by imposing the preserva- 
tion of the specific energy balance between the two major 
halo components during the accretion phase, but use in- 
stead the Lx — Tx relationship to adjust simultaneously 
the values of this parameter and of the excess energy, 
we would have inferred, precisely, a best-fit value of 1.2 
for the universal polytropic index. This argues in favor 
of the consistency of the X-ray data with the inside-out 
growth of the structure of galaxy systems within a uni- 
versal preheating scenario in which AT is about half a 
keV per particle. 

On the other hand, the fact that our model shows the 
closest agreement to the observations across several or- 
ders of magnitude in X-ray luminosity for a value of the 
preheating energy that compares favorably to constraints 
found in the literature is very assuring and can be viewed 
as a further endorsement of the consistency of our ap- 
proach. Thus, our best-fit value of 0.55 keV particle -1 
for this parameter agrees very well with the energy bud- 
gets inferred from recent analytical and numerical IHM 
models (Voit et al. 2002; Finoguenov et al. 2003), as well 
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Fig. 3. — Left: Comparison between our best solution (solid curve) with the predictions for an isothermal (7 = 1; short-dashed curve) 
and nearly adiabatic (7 = 1.5; dot-short-dashed curve) gas endowed with the same amount of injected energy. Right: Effects of metallicity 
on the luminosity-temperature relation. The two lowest curves are from a 7 = 1.2, AE = 0.55 keV particle -1 model with zero (short- 
dashed line) and solar (solid line) gas metallicities. The two upper curves are the corresponding predictions resulting from a non-radiative 
(7, AE) = (1,0) model (same coding as before). The meaning of symbol types is the same as in Fig. 2. Note how the curves from the 
model with null energy injection clearly fail to match the observational data. 



as with the results of attempts at directly measuring the 
extra energy injected into the gas (Lloyd-Davies et al. 
2000). To our knowledge ours is among the universal pre- 
heating models with the lowest heat input requirements. 
In this respect, it is worth noting that excess energies of 
only half a keV per particle might not be out of the reach 
of the total energy released into (though not necessarily 
retained by) the intergalactic medium by galaxy winds 
in starburst galaxies or by active galactic nuclei (AGN), 
to name only a few of the most popular heating sources, 
which comes to several keV per particle when averaged 
over all baryons in the universe (e.g., Cavaliere, Lapi, & 
Menci 2002). 

3.3. Sensitivity to Model Parameters 

While specifying AE completely determines the gas 
profiles, it is reasonable to wonder how our predictions 
depend on the exact values adopted for this parameter, 
for the fixed-by-theoretical-arguments polytropic index, 
as well as for the other parameters of the model that have 
been prefixed for definiteness. In the following we com- 
ment on this issue, focusing on the Lx — Tx relationship. 

As was already noted on Section 2.2.2, the effects of 
the energy injected into the gas at preheating become 
progressively noticeable with decreasing halo mass (tem- 
perature), reflecting the mounting contribution of AE 
to the total energy of this component. In particular, 
we find that the behavior of the predicted luminosity- 
temperature relationship below Tx ~ 2 keV becomes 
extremely sensitive to the exact value adopted for this 
degree of freedom, to the point that in the adjust- 
ment process the X-ray data allows one to easily dis- 
criminate between variations of AE larger than only 
0.05 keV particle" 1 in spite of the fact that, at these 
scales, the dispersion of the observations is substantial. 

While AE sets the normalization and convexity of 



the luminosity-temperature relationship at the low-mass 
end, the amplitude of this correlation at cluster scales 
depends essentially on the value of the polytropic in- 
dex. By varying 7 (at the cost of lessening the physi- 
cal consistency of the predictions), we have found that 
relatively small changes of one tenth around its pre- 
ferred value of 1.2 produce theoretical Lx ~ Tx curves 
that differ only marginally from the best solution. The 
results of Section 3.1 demonstrate, however, that this 
mild degeneracy vanishes completely when the preserva- 
tion of the balance between the total specific energies 
of the gas and dark-matter components after halo for- 
mation for the entire range of scales we are consider- 
ing is taken into account. In contrast, as the left panel 
of Figure 3 illustrates, when the extremal values of the 
polytropic index 7=1 (isothermal IHM) and 7 = 1.5 
(nearly adiabatic IHM) are combined with our best-fit 
value for the injected energy, the model predictions for 
the adopted ACDM cosmogony clearly fail to match the 
X-ray data. It is worth noticing, however, that our so- 
lutions for Tx < 2 keV are relatively insensitive to the 
value of 7, as far as it is larger than ^ 1.15. This would 
explain why injection models adopting 7 = 5/3 for the 
gas are successful in describing the low-luminosity group 
data in spite of the lack of observational support for an 
adiabatic gas distribution at these scales. Moreover, the 
fact that the isothermal solution runs essentially paral- 
lel to our best model across most of the cluster regime, 
implies that universal preheating models relying on the 
isothermality of the IHM can be matched up with the 
data if they are endowed with some additional freedom 
that allows one to renormalize, directly or indirectly, the 
X-ray luminosity (temperature). 

Other factors that can influence our model predictions 
are those related to the formation history of DM halos, 
such as the underlying cosmogony, specified by the val- 
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ues assumed for the cosmological parameters, and the 
adopted mass-concentration relation — we have verified 
that the predicted gas properties for our best model re- 
main essentially independent of the adopted halo forma- 
tion epoch even when an explicit dependence on it is 
introduced by solving equation (14) at Zf OT . With regard 
to the dark and baryonic matter contributions to the 
density parameter, any variation in f2b and/or f2 m that 
lowers the cosmic baryon ratio would lower the amount 
of baryons in the IHM and thereby the halo luminosities 
(and, of course, the hot gas mass fraction; see § 4.1). 
At the same time, the normalization of the perturbation 
spectrum, as, has repercussions on the convexity of the 
Lx — Tx relationship at cluster scales, which increases 
when the value of this parameter is lowered and vice 
versa. We also want to emphasize the fact — already 
stressed by Balogh et al. (1999) — that, in models re- 
lying on the hydrostatic equilibrium of gas in the dark 
matter potential, the ratio of model to observed luminos- 
ity is h dependent, so the best-fit solution varies with the 
assumed value of the Hubble constant. We have nonethe- 
less checked that if we take h — 1/2 then the preferred 
value of AE varies by less than 10%, provided we al- 
low for variations on the other cosmological parameters 
within the WMAP uncertainties. Last but not least, we 
also have verified that the adoption in our calculations 
of the mass-concentration prescription by Bullock et al. 
(2001) produces essentially identical results. 

The non-statistical noise of the Lx — Tx and any other 
scaling relationships involving bulk properties of the IHM 
may be increased by the intrinsic scatter of the mass- 
concentration relationship (e.g., Afshordi & Cen 2002; 
see also Fig. 1 in Wu et al. 2000), as well as by devi- 
ations of real systems from the ideal conditions of the 
modeling (e.g., non- negligible ellipticity, incomplete dy- 
namical relaxation, or variations in the efficiency of en- 
ergy injection). While these are effects that are not con- 
sidered in the present study, the gas metallicity is a fac- 
tor that may introduce a substantial amount of scatter 
in the L x — Tx relationship, especially at group scales, 
which is very straightforward to investigate within the 
context of our model. It is well established (e.g., Ren- 
zini 1997; Allen & Fabian 1998; Fukazawa et al. 1998; 
Ettori et al. 2001; Tozzi et al. 2003) that rich clusters 
exhibit almost constant average metallicity of <~ 0.3 Z Q 
with a reduced scatter, while cooler systems with T < 1 
keV show a wide range of values going from negligible 
metal contents to solar abundances — we do not consider 
here the existence of abundance gradients in cluster and 
groups of galaxies; see, e.g., De Grandi et al. (2004) and 
references therein. The fact that metal lines begin to 
provide a significant contribution to the emissivity only 
at temperatures below a few keV has been frequently 
exploited in investigations of the X-ray properties of ha- 
los by computing the cooling function under the simpler 
conditions that describe the physics of high temperature 
plasmas (but see Dave et al. 2002). However, when one 
seeks to extend the predictions to the group scales, the 
effects of metallicity can no longer be overlooked. As we 
show in Figure 3 (right), the contribution of metal lines 
to cooling starts to be noticeable, just like it happens 
with the contribution of AE to the gas energy, when 
the IHM temperature drops under 2-3 keV, to the point 
that, independently of the amount of non-gravitational 



heating, the bolometric X-ray luminosities at Tx ^0.5 
keV inferred from the model solutions that use a solar- 
mctallicity cooling function become about three times 
higher than those in which primeval abundances have 
been adopted. (The comparison between the two sets of 
preheated and non-preheated solutions drawn in the fig- 
ure also illustrates the strong downward bending in the 
luminosity-temperature relationship that results from an 
energy injection of 0.55 keV particle -1 , which for 0.5 keV 
halos implies a reduction of more than two orders of 
magnitude in Lx for 7 = 1.2.) Also note how within 
the isothermal models (7 = 1) with no energy injection 
(AE = 0), the upper curve, which represents the solution 
for a gas with a solar metal abundance, veers away from 
the roughly self-similar scaling of the zero-metallicity so- 
lution, showing a shallower, rippled slope that reflects the 
rises and falls of the cooling function arising from differ- 
ent metal lines (see, e.g., Figures 8 and 18 in Sutherland 
& Dopita 1993). 

Lastly, we would like to comment on the fact that some 
analytical models also resort to surface brightness bias, 
i.e., to measure temperatures and luminosities above an 
arbitrarily fixed surface brightness level, to explain both 
the offset toward lower luminosities and the increase of 
the dispersion of the luminosity-temperature relation for 
the less massive halos (e.g, Voit et al. 2002). Given the 
outward-decreasing equilibrium profiles of the gas, any 
reduction in the outer halo boundary has the effect of 
lowering Lx and raising Tx, producing the same trend 
seen in the data. Although this effect has not been ac- 
counted for in our modeling to avoid introducing an ad- 
ditional, poorly constrained, degree of freedom, we point 
out that raising the surface brightness threshold would 
have reduced further the amount of preheating energy 
necessary to reproduce the observations. 

Given that the (relatively important) impact of possi- 
ble variations in metallicity is noticeable for a range of 
system temperatures in which X-ray data — other than 
the X-ray luminosities and temperatures — are scarce 
and affected by large uncertainties, the discussion of the 
next section focuses only on the model predictions in- 
ferred for halos whose IHM has a metallicity fixed to the 
conventional value of 0.3 Zq. 

4. RESULTS 

4.1. The mass-temperature relationship 

Plotted in Figure 4 is a comparison between the mass- 
temperature relationship resulting from our preferred 
model and data from both the enlarged HIFLUGCS sam- 
ple (Reiprich & Bohringer 2002) and the Birmingham- 
CfA cluster scaling project (Sanderson et al. 2003). The 
temperature axis shows, as before, the emission- weighted 
temperature within the virial radius, Tx, while to fa- 
cilitate comparison with observations, the halo mass is 
now represented by M 50 o, the total mass inside a scaled 
radius of c 50 o within which the mean matter density 
is 500p cr it- As Sanderson et al. (2003) actually quote 
emission-weighted temperatures within a scaled radius 
of 0.3c2oo their values have been conveniently rescaled by 
means of the relationship Tx = Tx(0.3c2oo)/l-25, which 
we have calibrated from the comparison between their 
estimates of the emission-weighted temperatures for the 
subsample of Lloyd-Davies et al. (2000) 's systems and 
the measurements published in this previous study. 
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Fig. 4. — Mass-temperature relation versus observations. Differ- 
ent symbols are used to identify different datasets. The lower solid 
line shows the predictions corresponding to the best-fit model in 
Fig. 2. The two upper curves illustrate that the correlation between 
halo mass and gas temperature inferred in the AE = simulations 
of Evrard et al. (1996) for a IHM obeying an isothermal (7 = 1) 
/3-model with an iron abundance of 0.3 solar (long-dashed line) 
is closely matched by our corresponding model prediction (short- 
dashed line). 



The visual inspection of Figure 4 shows that our pre- 
dicted mass-temperature relationship has the right nor- 
malization and that the constant logarithmic slope of 
~ 1.7 of the model solution above ~ 1 keV is fully compat- 
ible with the trend shown by the observations. This good 
match between the observed and theoretical M500 — Tx 
relations is very compelling, especially considering that 
there is no freedom in our prediction. Our results are also 
in excellent agreement with previous findings of a loga- 
rithmic slope of 1.79 ± 0.14 by Nevalainen et al. (2000), 
by using overdensities of 1000 and spectroscopically de- 
rived temperatures, and of 1.78 ± 0.09 by Finoguenov 
et al. (2001), derived from spatially resolved X-ray ob- 
servations. In contrast, the slope of 1.84 ± 0.06 found 
by Sanderson et al. (2003) for their sample (using M 2 oo) 
differs by 2.3a from our value. However, a major dif- 
ference of this latter sample, compared to those used in 
the studies previously mentioned, is the inclusion of a 
substantial number of systems with temperatures below 
1 keV, including two measurements on early-type galax- 
ies (the lower-leftmost data point in Fig. 4 represents 
the most massive of them, NGC 6482), which tend to 
increase the slope of a global log- log linear fit. On the 
other hand, the inferred uniform rising of the mass over 
most of the temperature range does not support claims of 
a mass-temperature relationship that is convex (Sander- 
son et al. 2003; but see our comments below on the to- 
tal gas mass-temperature relationship) or has a break at 
about 3-4 keV (Nevalainen, Markevitch, & Forman 2000; 



Finoguenov et al. 2001). At this point, we would like to 
stress that any conclusions on the actual value of the 
slope of the mass-temperature relationship and its pos- 
sible convexity should be regarded with caution, mostly 
because the results being debated are not only sometimes 
based on ill-defined, incomplete datasets, but also due 
to the different prescriptions, extrapolations, and other 
inherent difficulties involved in the calculation of inte- 
grated halo properties, particularly the virial mass. The 
definitive answer will have to await the forthcoming avail- 
ability of larger samples of high-quality data with well- 
defined selection procedures which can be accounted for 
during the fitting process. What is clear from our anal- 
ysis is that the predictions of our universal preheating 
model compare very favorably to current observations. 

For low-temperature systems (Tx < 1 keV) the pre- 
dicted M500 — Tx relationship initiates a break away from 
the general M 5no oc T£ 1-7 trend, barely perceptible be- 
cause of the cut off in mass (and temperature) imposed 
by the condition that the IHM must be bound for excess 
energies of 0.55 keV particle -1 . Had we not taken into ac- 
count this restriction the theoretical mass-temperature 
relation would had then shown a progressively marked 
decline, leading to an emission-weighted temperature of 
the gas almost independent on the halo mass, much as 
the universal preheating model by Balogh et al. (1999) 
and Babul et al. (2002) predicts it should happen for 
a range of masses M 500 around 10 13 M Q . Most interest- 
ingly, this radical variation of the correlation trend would 
show a very good correspondence in the diagram with the 
locus occupied by the measurements of the two early-type 
galaxies included in the Birmingham-CfA cluster sample. 
This coincidence is, nonetheless, surprising since, for the 
best-fit value of the preheating energy, our model gives 
rise to galaxy-sized halos that cannot retain their hot 
gas. 

Concerning the pure preheating model by Balogh, 
Babul, and collaborators, we want to stress the funda- 
mental difference existing in the explanation put forward 
for the origin of the observed dispersion in the scaling 
relationships, which these authors attribute primarily to 
the distribution of halo formation times. This conclu- 
sion appears to be a consequence of their assumption 
that the accretion of the gas on the low-mass halos is 
specified by the adiabatic Bondi accretion rate, which 
leads to a strong dependence of the gas density normal- 
ization on Zf or that propagates, throughout the Jeans 
equation, to the bulk properties of this component. In 
contrast, the properties of the X-ray gas in our best 
model are independent of the distribution of halo for- 
mation redshifts. Therefore, taking into account both 
that, on cluster scales, the other intrinsic dispersion fac- 
tors that might affect the group data are not an issue 
and that the level of scatter of the measurements de- 
picted in Figure 4 is pretty uniform across the range of 
observed masses, we attribute the observed scatter in the 
-^500 — Tx relationship a genuinely statistical nature. 

We have also included in Figure 4 the self-similar re- 
lation derived from the hydrodynamical simulations by 
Evrard et al. (1996), rescaled to our value of h (long- 
dashed curve). Both the higher mass normalization and 
the canonical 1.5 slope are faithfully described by our 
model prediction corresponding to an isothermal gas 
without preheating (short-dashed curve). This result, 
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Fig. 5. — Gas mass-temperature relation versus observations. 
Solid circles and error bars represent gas mass measurements within 
a scaled radius of C500 by Mohr et al. (1999). The lower solid line 
shows the prediction corresponding to the best-fit model in Fig. 2. 
The upper short-dashed line shows the prediction resulting from a 
model with 7 = 1 and AE = 0. 



Fig. 6. — Relation between the integrated gas mass fraction and 
gas temperature. The data points are gas mass fraction measure- 
ments within a scaled radius of C500 from Mohr et al. (1999). The 
solid curve is the prediction corresponding to the best-fit model in 
Fig. 2. 



apart from reinforcing confidence in our model, demon- 
strates that essentially all the steepening, as well as the 
downward shift in the normalization, are attributable to 
the effects of preheating. 

The impact of preheating on total gas mass obser- 
vations can be gauged in Figure 5, where we depict 
emission-weighted temperature versus gas mass -Mx,5oo 
inside the scaled radius c 50 o- Again, the trend exhib- 
ited by the data, now coming from Mohr, Mathiesen, & 
Evrard (1999), is dutifully reproduced by our relation- 
ship. Another remarkable result is that we predict a 
convex correlation in which the slope smoothly increases 
with decreasing temperature as in the purely analyti- 
cal model by Dos Santos & Dore (2002). The compar- 
ison of this correlation with the result for an isother- 
mal model in which preheating has been switched off 
(short-dashed curve), shows that both predictions follow 
parallel trends for systems with Tx > 5 keV (the nor- 
malization of the isothermal model always being higher). 
For halos this hot, the inferred typical logarithmic slope 
is ~ 1.7, close to the self-similar value but nevertheless 
steeper, and an almost perfect match to the best-fit slope 
of 1.71 ± 0.13 found by Vikhlinin et al. (1999) for regu- 
lar galaxy clusters. Below this temperature, as the role 
of preheating becomes progressively important, the gas 
mass-temperature relationship of our best model tends to 
steepen promptly, adopting power-law indexes between 
2.3-2.5. Not surprisingly, the best-fit linear log- log cor- 
relation determined by Mohr et al. (1999) to their data 
produced a slope of 1.98 ± 0.18, intermediate between 
these two behaviors. 



The fact that our (■y,AE) = (1,0) solution exhibits a 
nearly constant slope (there is only a hint of convexity) 
somewhat steeper than the self-similar value of 1.5 does 
not come as a surprise if one considers that the behav- 
ior of the Mx.500 — Tx relationship reflects those of the 
ratios Tx/T viv and M x /M. While the identity between 
the effective gas and dark matter virial temperatures and 
the proportionality of the density distributions of these 
two main components for all radii are common assump- 
tions in the modeling of X-ray halos, our results illustrate 
that they are only convenient approximations. Thus, 
within our scheme, even an isothermal gas with null non- 
gravitational heating (otherwise the gas temperature is 
pushed well above T v j r ), leads to values of the ratio be- 
tween the X-ray luminosity-weighted and virial tempera- 
tures that, in general, border on, but are not equal to the 
unity (indeed, in all the solutions investigated this ratio 
shows a monotonic increase with increasing halo concen- 
tration). On the other hand, the hotter systems tend to 
be more gas-rich than the cooler ones, a tendency that 
preheating accentuates by preventing gas from reaching 
a high density in the central regions (Bialek et al. 2001; 
McCarthy et al. 2002; see also Sec. 4.3.1). 

Finally, we have compared in Figure 6 the temperature 
dependence of the gas mass fraction as predicted by our 
preferred model — calculated by simply dividing the pre- 
viously inferred masses Mx.500 and M500 as a function 
of the emission-weighted temperature — with data from 
Mohr et al. (1999). Looking at this plot one sees that, 
regardless of the considerable amount of scatter in the 
measurements, the prediction corresponding to our pre- 
ferred model again explains the observations remarkably 
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Fig. 7. — Left: Inner halo entropy at one tenth of r2oo plotted against the system mean temperature within 0.3r2oo- Data points are 
measurements from Ponman et al. (2003). The right- most solution (short-dashed line) indicates the effect of switching off preheating (i.e., 
taking AE = 0). The entropy floor of 126keVcm 2 from Lloyd-Davies et al. (2000) is indicated by a thin horizontal dotted line. Right: 
Scaled outer halo entropy at the overdensity of 500 vs. total mass measured within the same overdensity. Data points show measurements 
compiled by Finogucnov ct al. (2002). The near horizontal short-dashed curve shows the solution when AE = 0. In both panels the curves 
with the largest amplitude (solid lines) are the predictions corresponding to the best-fit model in Fig. 2. 



well. One can also observe that our model prognosti- 
cates a progressive steepening of the slope, from cluster 
to group scales, which we cannot validate because of the 
insufficient dynamical range of the data available (but 
note the apparent consistency of our prediction with the 
trend delineated by the data of Sanderson et al. 2003 
in their Fig. 5). We find that the total gas mass frac- 
tion within C500, -Px,500) levels off to values within about 
±10% of the cosmic value for emission weighted tem- 
peratures above 3-4 keV (i.e. for virial masses above 
(3-4)xlO 14 M0), while it becomes directly proportional 
to Tx for less massive systems, which, in turn, implies 
that -Fx, 500 oc M®qq in the hot group regime. This is a 
mass dependence somewhat weaker than the self-similar 
behavior Fx oc M predicted by the models of Balogh 
et al. (1999) and Dos Santos &' Dore (2002). Note also 
that our results are consistent with the well-known fact 
that pressure forces arising from a high entropy floor (see 
next section) can efficiently shut off the gas supply to 
the halos, reducing the amount of baryons that end up 
gravitationally bound (Cavaliere et al. 1998; Muanwong 
et al. 2002; Oh & Haiman 2003) and resulting therefore 
in smaller total baryon fractions for halos with shallow 
potential wells. 

4.2. Scaling Relations for the Entropy 

A particularly interesting source of information about 
deviations from self-similarity is the gas entropy, as it 
is a record of the thermodynamic history of the diffuse 
baryons in clusters and groups of galaxies. In X-ray 
studies this property is customarily represented by the 
entropy of the gas electrons, which we infer from the ex- 
pression 



with ((T, Z) the function that renormalizes the number 
density of gas particles to the number density of elec- 
trons 5 . In the last few years, a substantial amount of 
observations on the scaling properties for the entropy 
on both the inner and outer regions of galaxy systems 
have become available. As noted earlier in the Intro- 
duction, such observations have served to establish that 
X-ray clusters and groups show entropy excesses with re- 
spect to the expectations from pure gravitational shock- 
heating (e.g., Ponman et al. 1999; Lloyd-Davies et al. 
2000) that are not restricted only to their central regions 
(Finoguenov et al. 2002; Ponman et al. 2003). 

In order to confront our predictions with data from the 
observational studies just mentioned, we have calculated 
the central gas entropy (from eq. [19]) at a scaled radius 
of 0.1c 2 oo against mean gas temperature within 0.3c2oo, 
while to estimate the entropy level on the halo outskirts 
we have chosen the adiabat at C500 scaled by M 500 ' 13 , the 
total mass within this radial distance expressed in units 
of 10 13 M Q , as a function of the total system mass M 50 o- 

Data and model predictions are compared in Figure 7. 
For the inner entropy (left panel), our results describe, 
reasonably well, the data points from Ponman et al. 
(2003) given their large observational uncertainties. As 
shown in the plot, our prediction deviates progressively 
from the K oc T scaling followed by the high tempera- 
ture systems, reaches a minimum value between 1 and 
2 keV, and rises again for the coldest objects. Remark- 
ably, this behavior is relatively consistent with the trend 
described by the data in the pioneering studies of the 
entropy by Ponman et al. (1999) and Lloyd-Davies et al. 
(2000), which predicted an entropy floor at group scales 
with a typical value of 126kcVcm 2 (for our adopted h) 



5 As in our model 7^5/3, the ratio T/n^ 3 depends, for a given 
M, on the haloccntric distance. 
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not too different from the minimum value reached by our 
model, although in their sparse data there was no sign 
of a recovery of the entropy at the low-temperature end. 
However, the measurements in the much larger sample 
gathered by Ponman ct al. (2003) appear rather to delin- 
eate an 'entropy ramp', meaning that the departure from 
the self-similar scaling is more accurately represented by 
a gradual monotonic deviation characterized by a loga- 
rithmic slope of about 2/3. Quite interestingly, though, 
further data coming from elliptical galaxies (O'Sullivan, 
Ponman, & Collins 2003) — which we recall cannot be 
consistently explained by our modeling — show some ob- 
jects with very high central entropies. Given the good 
agreement of the predicted gas temperatures with obser- 
vations (see Fig. 4), this possible mismatch of the central 
entropy for the coldest systems may be attributed to the 
fact that their associated central gas densities arc typi- 
cally lower than actual measurements (compare, for in- 
stance, our results in Sec. 4.3.1 with Fig. 8 of Sanderson 
& Ponman 2003). 

In Figure 7 (right) the predictions of our preferred 
model for the scaling of the 'reduced' outer entropy, 

^x(c5oo)/-^5oo i3> against system mass are compared 
with the data from Finoguenov et al. (2002). As seen 
from the figure, in this case the mean trend drawn by 
the full set of data points is followed with impressive fi- 
delity by the theoretical curve. Nevertheless, our results 
introduce two major modifications with respect to the in- 
terpretation given by Finoguenov and coworkers. First, 
we find that excess entropy should be present across the 
full mass range of galaxy systems and not only in those 
below a certain mass threshold (in line with the inter- 
pretation given in Ponman et al. 2003). And second, 
the scaled entropy of the model shows a monotonic in- 
crease with decreasing halo mass, thereby implying the 
absence of the 'entropy ceiling', i.e., of the upper limit 
on the preheating, that was claimed to be present in 
the ASCA observations. As in Finoguenov et al. (2002), 
we also have included in the plot the entropy resulting 
from purely gravitational heating (accomplished by set- 
ting AE to zero). Under these conditions, we predict 
a nearly constant scaled entropy of <~ lOOkeVcm 2 - 
consistent with the expectation that under self-similarity 

conditions scaling by M 500 ' should renormalize the en- 
tropy to a value independent of system temperature — 
that clearly fails to explain the observations. 

4.3. Equilibrium Configuration of the Gas Distribution 

We now present our model predictions for the inter- 
nal structure of the IHM. We first focus on the spatially 
resolved properties of the X-ray-emitting gas and then 
confront the gradients of the surface brightness and pro- 
jected temperature profiles with X-ray data. 

4.3.1. Intrinsic Profiles 

In Figure 8 we plot the three-dimensional radial profiles 
of the scaled X-ray gas temperature, mass density, and 
entropy, as well as of the integrated gas mass fraction, 
for a set of representative halo masses at z = 0. Radial 
distances are given in units of r v i r so as to facilitate the 
comparison among profiles of different masses and with 
the results of previous works. Quite remarkably, the gas 
temperature and density profiles (and hence the entropy 
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Fig. 8. — Predicted three-dimensional radial profiles of the X- 
ray gas temperature, density, entropy, and integrated gas mass 
fraction for five present-day halos of different total masses. The 
three top panels show profiles scaled to their central values, while 
the cumulative baryon fraction profile is given in units of the cosmic 
value. Radial distances are normalized to r v ; r . The straight short- 
dashed line included in the third panel is used to indicate the slope 
of 1.1 expected for the scaled entropy from pure shock heating. 



profiles) for halo masses above 6.3 x 10 13 M Q look very 
similar when scaled to their central values, although we 
recall that, in profiles of the NFW-like form, the radial 
coordinate must be expressed in r s units in order to make 
it truly independent of the halo mass (see, for instance, 
the model eqs. in § 2.1 and 2.2). Only the less massive 
(coldest) systems tend to have slightly shallower profiles 
as a result of the increasing influence of preheating. In 
contrast, the bottom panel of Figure 8 shows that the 
variation of the predicted gas fraction with radius, which 
with the exception of the outer regions (x > 0.4c) of the 
most massive systems (M > 6.3 x 10 14 M Q ) decreases 
monotonically toward the halo centers, depends strongly 
on the halo mass (temperature), mirroring the trend seen 
in Figure 6 for the global values: the more massive the 
halos, the higher the profile amplitude. This behavior 
agrees quite nicely with the most recent observational de- 
terminations of this property by Sanderson et al. (2003), 
illustrating that the distribution of the IHM does not 
follow that of the dark matter, but is significantly less 
concentrated as was expected (see below). 

The good agreement with observations is specially ex- 
citing for the entropy profiles. On one hand, none of our 
profiles, even for the smallest halos, is fully iscntropic, 
thus invalidating one of the major criticisms against uni- 
versal preheating gas models consisting in the production 
of flat entropy gradients for group-sized halos. On the 
other hand, our model reproduces quite well, for all but 
the less massive systems, the outer radii K(r) oc r 11 be- 
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Fig. 9. — Comparison of the predicted X-ray gas density profiles 
resulting from our preferred model (thick curves) and those result- 
ing from a model with 7 = 1 and AE = (thin curves). Each pair 
of equal-type curves corresponds to a halo of a given current total 
mass. Radial distance is in units of r v i r . 



havior observed by Ponman et al. (2003) and predicted 
in the gravitational shock-dominated regime of entropy 
production by both theoretical prescriptions based on 
spherical accretion within a NFW dark matter halo (Eke 
et al. 1998; Tozzi & Norman 2001; Babul et al. 2002) 
and the latest state-of-the-art cosmological gasdynami- 
cal simulations (Borgani et al. 2004). This result indi- 
cates that it is possible to reconcile the large entropy 
excesses seen at large radii, which suggest that entropy 
profiles are dominated throughout by the effects of non- 
gravitational preheating, with the fact that their typi- 
cal outer slopes scatter about the value predicted from 
shock heating, without requiring any specific fine tunning 
of the model. Interestingly enough, the simulations by 
Finoguenov et al. (2003) show that heating with an equal 
amount of energy per particle produces entropy profiles 
that are similar to those arising from accretion shocks of 
fixed strength. 

A final aspect concerning the entropy profiles that we 
wish to point out is that, if we apply the same scalings 
as in Finoguenov et al. (2002), our universal preheating 
model reproduces, very satisfactorily, the behavior and 
amplitude at all radii of the deprojected profiles mea- 
sured by these authors (see also Ponman et al. 2003). In 
particular, we are capable of replicating the strongest rise 
with radius of the richest clusters and both the flatter and 
higher entropy levels, especially at smaller radii, found 
for groups. As we have not included radiative cooling in 
our treatment, this suggests that the possible effects of 
this latter mechanism might not be excessively impor- 
tant, even in the central regions of massive clusters. 

The three-dimensional radial profiles of the gas den- 



sity have been also redrawn in Figure 9 to provide an 
example of the impact of preheating on the gas distri- 
bution. As in Figure 8, we plot the profiles, now ex- 
pressed in physical units to allow a direct comparison 
with observations, corresponding to five halos with dif- 
ferent masses ranging from 2 x 10 13 to 2 x 10 15 M . Two 
sets of curves are shown: one corresponding to our best 
model and the other to the canonical (7, A£) = (1,0) 
'non-radiative' conditions. It can be seen that preheat- 
ing manifests itself as a substantial reduction of the hot 
gas density in the central halo regions. In agreement 
with expectations, this reduction is strongest for the low- 
mass systems. Clearly, as the halo mass is lowered, the 
gas becomes less concentrated and exhibits a bigger core. 
Also evident is the fact that the gas distribution in the 
outer halo regions for systems with M > 6.3 x 10 13 M Q 
is more extended (i.e., the density at a given radius is 
raised) than that corresponding to the (7, AE) = (1,0) 
prediction. In contrast, the effects of preheating lead 
to a substantial increase in the gas temperature and en- 
tropy toward the halo centers (as similarly reported by 
McCarthy et al. 2002 and Babul et al. 2002), not quite 
unlike the predictions by galaxy formation-regulated gas 
evolution models (e.g., Wu & Xue 2002). 

4.3.2. Surface Brightness Profiles 

The observed X-ray surface brightness at a projected 
scaled (by r s ) radius X from the center of a spherical 
halo is given simply by the Abel integral 6 

S x (X) = 2r s e Av (x)di, (20) 

Jo 

where I = \J x 2 — X 2 is the distance along the line-of- 
sight and cav{x) is the X-ray volume emissivity within 
the energy band Av of the observations. It turns out 
that our predicted X-ray surface brightness profiles in 
the ROSAT broad (0.1-2.4 keV) band are a reasonable 
match to the conventional /3-model, Sx(R)/Sx{0) = 
[l + ( J R/r c ) 2 ]" 3 ' 3+1 / 2 , within c 500 . Beyond this radius, the 
projected gas emission initiates a rapid downfall due to 
the finite outer boundary of the integral (20). Thus, for 
the purposes of comparison with the existing X-ray imag- 
ing measurements, the resulting Sx profiles of present- 
day halos have been matched by /3-models over the radial 
interval O.OIC500 < X < c 50 o, which corresponds approx- 
imately to extending the fit out to half the virial radius, 
and the best-fit values of the slope (3 and core radius r c 
have been compared with their observational determina- 
tions. We note that the portion of the profile being fitted 
is consistent with the fact that measurements of the X- 
ray surface brightness of clusters and groups generally 
embrace only a fraction of their total radius. 

Figure 10 shows how the predicted dependence of our 
best-fit /3-model parameters on Tx compares to data ex- 
tracted from the same catalogs used in the investiga- 
tion of the mass-temperature relation, with the emission- 
weighted temperatures quoted by Sanderson et al. (2003) 
conveniently rescaled (see § 4.1). Although the large dis- 
persion in each parameter at a given temperature does 

6 While strictly the integral along the line-of-sight (20) should 
extend to infinity, the applicability of the Jeans equation (8) be- 
yond r v ; r becomes questionable. 
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Fig. 10. — Temperature dependence of the slope parameter (top) 
and core radius (bottom) inferred from f3- model fits to the X-ray 
surface brightness profiles. Different datasets are identified with 
different symbols. Curves show the results of fitting within the 
radial range O.Olcsoo < X < C500 the Sx in the 0.1-2.4 keV band 
predicted by our preferred model. 



little to constrain the model predictions, our theoreti- 
cal curves fare, once more, quite satisfactorily, as they 
appear to roughly track the mean trends that would be 
inferred from 'chi-by-eye' fits to the data points. Quite 
remarkably, the prediction for the slope parameter sup- 
ports the frequently advocated case for flatter surface 
brightness profiles for groups (e.g., Ponman et al. 1999; 
Helsdon & Ponman 2000; Mulchaey et al. 2003; Osmond 
& Ponman 2004), while for high-mass systems it typically 
approaches the self-similar value of 2 /3 (Jones & Forman 
1984). On the other hand, the predicted increase of the 
central gas density concentration with the ~3/4 power 
of the gas temperature is also consistent with the data, 
which suggests further a possible sharp drop of r c for sys- 
tems with X-ray temperatures < 1 keV. Regarding this 
latter possibility, we point out that, as shown by Voit et 
al. (2002), observational systematics such as the range of 
radii in the fit and surface brightness bias can affect mea- 
surements for low temperature systems. In the former 
case, we have verified that fitting over the radial range 
0.01c 50 o < X < 0.5c 50 o essentially involves an overall re- 
duction of the amplitude of the predictions, while raising 
the surface brightness threshold, which makes the best-fit 




Fig. 11. — Predicted projected emission-weighted X-ray gas tem- 
perature profiles resulting from our best model, scaled by Tx,5oo- 
Thick curves represent five halos of different total masses. The 
dotted boxes enclose the ASCA temperature profiles observed by 
Markevitch et al. (1998): the outer one shows the 90%-error band 
of all their temperature profiles, while the inner box approximates 
the scatter of the corresponding best-fit polytropic models. Aver- 
age temperatures of the observed clusters range from 3 to 10 keV, 
which in our model correspond to total halo masses ranging from 
~ 2 X 10 14 to 2 X 10 15 M Q . To facilitate comparison with data 
projected radial distance is in units of riso- 



values of (} and r c both decline faster with decreasing Tx 
has, however, the undesirable consequence of excluding 
from detection the coldest groups whose distinct behav- 
ior one is trying to describe. 

4.3.3. Projected Temperature Profiles 

As a last check of the consistency of our model, 
we have calculated the two-dimensional projection of 
the emissivity-weighted temperature profile of the halos 
(eq. [17]), which is defined as 



Tx(X) 



V 2 P 

v c Jo 



- A 2 



fl(T* Z)T*(x)e{x)dl 



V7: 



-X' 1 



e(x)dl 



(21) 



Figure 11 illustrates the predictions by our preferred 
model for the projected emission-weighted temperature 
profiles of five halos with present-day masses equally 
spaced (in logarithmic scale) between 2 x 10 13 and 
2 x 1O 15 M0, which rise uninterrupted toward the halo 
centers. Interestingly, our results for masses above 
2 x 1O 14 M (i.e., T x > 3 keV) compare very favor- 
ably with the composite radial temperature profile by 
Markevitch et al. (1998). From an analysis of ASCA spa- 
tially resolved spectroscopic data for a set of nearby clus- 
ters with mean spectral temperatures above 3 keV, these 
authors found that the radial run of the temperature 
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for symmetric, cooling-flow corrected, clusters showed 
a nearly universal outward decline between ~ 0.1 and 
~ 0.6ciso corresponding to a polytropic index 1.24+Q.12 
(with 90% errors). We note that the model profiles have 
been scaled by Tx.sooj the average temperature within 
C500, in an attempt to represent more accurately the scal- 
ing by the average temperatures calculated over the ra- 
dial range covered by Markevitch et al.'s data. A differ- 
ent (similarly realistic) normalization would have only a 
marginal impact on the profiles and would not change 
our conclusions. 

Albeit in another work based on ASCA observations 
Finoguenov et al. (2001) reported a similarly outward de- 
creasing behavior, the results by Markevitch et al. (1998) 
have nonetheless been controversial, with some other 
ASCA papers claiming consistency with overall isothcr- 
mality (White 2000). Besides, there is the well-known 
analysis of Beppo-SAX data by De Grandi & Molendi 
(2002), who derived profiles similar to those of Marke- 
vitch et al. (1998) in the outer cluster regions, but not 
in the very centers, where they observed flat tempera- 
ture gradients for radii smaller than ~ 0.2ciso- Other 
Beppo-SAX studies of a more limited radial extent pro- 
duced temperature profiles that were generally flat or 
even rising slightly out to ~ 30% of the virial radius (Ir- 
win & Bregman 2000). Now, however, the latest results 
by Chandra appear to confirm quite nicely the steady de- 
crease of the gas temperature in the 0.1 — 0.6ciso radial 
range, although the central peak at x < O.lciso looks like 
an artifact of the cooling flow modeling, as ASCA could 
not resolve the centers well 7 (Markevitch 2004, private 
communication). Further support to an outward decreas- 
ing IHM temperature over a great deal of the radial range 
comes from the results of gasdynamic simulations both 
adiabatic (e.g., Ascasibar et al. 2003) and with cooling, 
star formation, and supernova feedback (Borgani et al. 
2004; Ettori et al. 2004). 

5. SUMMARY AND CONCLUSIONS 

We have investigated the global and structural X- 
ray properties of the IHM in nearby groups and clus- 
ters of galaxies within a flat ACDM cosmology using 
a simple, reasonably well specified analytic polytropic 
gas model relying on the inside-out growth of structure 
between major mergers. Among the different scenarios 
suggested in the literature as a means of explaining the 
well-established breakdown from the simplest self-similar 
predictions, we have focused on the one which assumes 
that the specific energy of the gas is increased by some 
arbitrary source of extragravitational heating before the 
virialization of galaxy systems. This additional injection 
of energy onto the intergalactic medium has been pre- 
sumed to be independent of the virial mass of the halo 
where the hot gas is incorporated. We summarize our 
main findings below. 

(i) A polytropic index 7 = 1.2, as many observations 
and numerical simulations suggest, offers the best consis- 
tency with the postulate that the specific energy balance 
between the two main halo components (dark matter and 

7 As our model, by design, ignores radiative cooling, we cannot 
account for the smooth decline of the temperature toward the cen- 
ter observed in the innermost regions (< O.lcigo) of some massive 
clusters with strong cooling flows (e.g., Allen, Schmidt, & Fabian 
2001). 



hot gas) is preserved when energy loses are negligible. 
This is the first time a theoretical model tackling the X- 
ray properties of relaxed galaxy systems — the condition 
that the preheated gas remains bound constraints our 
predictions to halos with present-day total virial masses 
above ~ 2 x 10 13 M — provides a possible physical origin 
for such a value of 7. A natural consequence of this hy- 
pothesis is that the IHM properties must be essentially 
independent of the system formation redshift. Variations 
on the metallicity of the hot gas could be a possible alter- 
native to claims that favor differences in the formation 
epoch as the primary source of the dispersion in the ob- 
served scaling relations at group (Tx < 2 keV) scales. 

(ii) With 7 fixed to 1.2, our hot X-ray gas model for 
group- and cluster-sized halos is capable of describing 
very accurately the observed mean trends of a set of rep- 
resentative scaling laws between bulk properties of this 
constituent for a level of energy injection of 0.55 keV 
per gas particle, low enough to be accessible to popular 
nongravitational heating sources such as supernovae and 
AGN. In particular, we infer mass-temperature relation- 
ships with slopes steeper than the self-similar prediction. 
For the M500 — Tx relationship we find a nearly constant 
logarithmic slope of ~ 1.7 across most of the temperature 
range that initiates a slight offset toward lower masses for 
Tx < 1 keV pointing to the locus occupied by early-type 
galaxy data. On the other hand, the Mx,5oo — Tx corre- 
lation exhibits a slope similar to the M500 — Tx relation 
for the hottest clusters that progressively steepens as the 
scale diminishes, suggesting that for group-sized systems 
the total gas fraction within the halo should grow steadily 
with mass approximately as M 6 . The model also allows 
us to confidently rule out purely isothermal (7 = 1) or 
adiabatic (7 = 5/3) descriptions of the X-ray data. 

(iii) The scaling properties of the gas entropy are sim- 
ilarly well predicted. According to the results of our 
investigation, a preheating stage involving a relatively 
modest amount of energy injection of about half a keV 
per particle, suffices to explain, simultaneously, the near 
lOOkeVcm 2 level shown by the central entropy of halos 
at group scales and the higher entropy level of ~ 400- 
500kcVcm 2 seen in the periphery of these galaxy sys- 
tems. 

(iv) Only the prognosticated behavior of the inner en- 
tropy for the coldest systems appears to be in conflict 
with observations. Specifically, in halos with emission- 
weighted temperatures below ~ 1 keV our model indi- 
cates that the central entropy should increase with de- 
creasing gas temperature, probably due to the substan- 
tial reduction in their central gas densities caused by pre- 
heating. In contrast, the preferred interpretation of the 
most recent observational data is that the central entropy 

obeys a non-self-similar power law form, Kx oc T x , 
at all scales. However, the systematic and statistical 
uncertainties in observations resulting from limitations 
in the instrumental response and the intrinsic faintness 
of the groups are sufficiently important that a strong 
discrepancy cannot be claimed. Alternatively, one may 
also consider the possibility that a universal preheating 
model docs not provide an accurate description of the 
IHM physics in the innermost regions of cold galaxy sys- 
tems. This would be the case if the epoch of preheating 
occurs relatively late in the history of the universe, say 
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at z ~ 2, so that most of the group cores, which tend 
to form early (zf or ~ 3-4) and settle later in the central 
regions of larger units, would be already in place when 
the energy boost takes place (calculations by Oh & Ben- 
son 2003 using the Press-Schechter theory show that at 
z = 2.2 about 50% of the cluster progenitors would not 
be affected by preheating). This implies that heating 
could have been less efficient in raising the entropy of 
less massive systems because of the high density of their 
IHM, whose thermodynamic history would be then dom- 
inated by the small amounts of entropy acquired through 
gravitational shocks. 

(v) The equilibrium gas profiles also agree quite nicely 
with recent observational results, suggesting that the as- 
sumed inside-out growth of structure after virialization 
offers a straightforward and sensible approximation to 
the evolution of the gas equilibrium profiles. Of par- 
ticular relevance are our findings that the IHM does 
not show, contrarily to other pure preheating models, 
a largely isentropic core and that the outer logarithmic 
slopes of the entropy profiles for all but the coldest sys- 
tems are quite similar to the 1.1 value expected from 
shock heating — without necessarily implying that this is 
the basic mechanism responsible for generating the ob- 
served excess entropy at large radii in groups and clusters 
of galaxies. 

(vi) We infer, in addition, outward declining tempera- 
ture profiles, much like most of the observational stud- 
ies report outside the cluster centers and in line with 
the trends found in recent hydrodynamical simulations. 
Claims that the hot gas distribution within cluster-sized 
dark halos is well represented by a polytropic model 
with 7 ^1.2 are reinforced further by the behavior 
predicted for the radial distribution of the projected 
emission-weighted temperature in relatively hot clusters 
(i.e., Tx > 3 keV). However, our model cannot account 
for either the nearly isothermal profiles or the large flat 
cores detected in certain recent observational analyses 
of temperature gradients in clusters. Although the fail- 
ure to reproduce these trends could reflect the need for 
incorporating some additional physics, such as a heat 
transport mechanism (e.g., thermal conduction) capable 
of bringing large regions of the IHM to the same temper- 



ature in a relatively short timescale, the fact is that the 
actual form of the profile that sets the norm for bound 
galaxy systems has yet to be settled. 

Our conclusion from these results is that, in general 
terms, our universal preheating model predictions are 
very succesful in fitting the observational data. The 
good match is more remarkable if one takes into account 
that it has been achieved without resorting to elaborated 
schemes of energy injection and the additional freedom of 
including the effects of observational biases in the mea- 
sured gas properties. Admittedly, processes deliberately 
excluded from our modeling, such as the condensation 
and removal of the lowest entropy gas from the IHM, 
must happen at some level because that is how galaxies 
and stars form within dark halos. Moreover, radiative 
cooling can help to reduce the energetic requirements for 
preheating while it might be necessary in order to explain 
the behavior of the temperature profiles in the innermost 
regions (r < 0.1r v j r ) of massive cooling-flow clusters re- 
vealed by new high-spatial resolution X-ray data. Nev- 
ertheless, the fact that the picture we are proposing can 
account, within current observational constraints, for the 
typical X-ray characteristics of nearby galaxy systems ar- 
gues in favor of it as a useful description of the evolution 
of the intergalactic medium properties in the densest re- 
gions of the universe. 

Yet the most important asset of the present work is per- 
haps that the X-ray properties of galaxy systems ranging 
from the coolest groups to the hottest clusters have been 
intimately linked to a successful analytic theory of the 
clustering history in the universe, thus providing an in- 
tertwined treatment for the formation and evolution of 
the structure of both nonbaryonic and baryonic matter 
that can be easily implemented in models of galaxy for- 
mation. 
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